arXiv:1505.04554vl [physics.space-ph] 18 May 2015 


JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1029/, 


Generation of Alfvenic Waves and Turbulence in Reconnection 
Jets 

Masahiro Hoshino , 1 Katsuaki Higashimori , 1 

Abstract. The magnetohydrodynamic linear stability with the localized bulk flow ori¬ 
ented parallel to the neutral sheet is investigated, by including the Hall effect and the 
guide magnetic field. We observe three different unstable modes: a “streaming tearing” 
mode at a slow flow speed, a “streaming sausage” mode at a medium flow speed, and 
a “streaming kink” mode at a fast flow speed. The streaming tearing and sausage modes 
have a standard tearing mode-like structure with symmetric density fluctuations in the 
neutral sheet, while the kink mode has an asymmetric fluctuation. The growth rate of 
the streaming tearing mode decreases with increasing magnetic Reynolds number, while 
the growth rates of the sausage and kink modes do not depend strongly on the Reynolds 
number. The sausage and kink modes can be unstable for not only super-Alfvenic flow 
but also sub-Alfvenic flow when the lobe density is low. The wavelengths of these un¬ 
stable modes are of the same order of magnitude as the thickness of the plasma sheet. 

Their maximum growth rates are higher than that of a standard tearing mode, and un¬ 
der a strong guide magnetic field, the growth rates of the sausage and kink modes are 
enhanced, while under a weak guide magnetic field, they are suppressed. For a thin plasma 
sheet with the Hall effect, the fluctuations of the streaming modes can exist over the plasma 
sheet. These unstable modes may be regarded as being one of the processes generating 
Alfvenic turbulence in the plasma sheet during magnetic reconnection. 


1. Introduction 

Alfvenic fluctuations and turbulence with orders of mag¬ 
nitude from several tens Re down to ion inertia lengths are 
often observed in the earth’s magnetotail [e.g., Russell, 1972; 
Hoshino et al., 1994; Bauer et al, 1995; Zelenyi et al, 2014]. 
The amplitude of these wave fluctuations in the plasma sheet 
can reach up to half the value of the lobe magnetic field, and 
these waves can carry a significant fraction of the plasma 
energy in the plasma sheet [e.g., Borovsky et al., 1997]. It 
has been discussed that turbulence with a large wave energy 
density plays an important role in the transport of mass and 
momentum in the magnetotail [e.g., Borovsky and Funsten, 
2003; Zimbardo et al, 2010]. Understanding the effects of 
this turbulence is believed to be a key factor in obtaining 
an overall understanding of the various plasma phenomena 
that occur both in space, and also in astrophysical plasmas 
[e.g., Birn et al, 2012; Lazarian et al, 2012]. 

To date, many researchers have discussed the importance 
of turbulence in the plasma sheet in the context of mag¬ 
netic reconnection. It has been argued that turbulence plays 
a crucial role in the dynamic evolution of magnetic recon¬ 
nection and that the magnetic energy dissipation rate can 
be enhanced by turbulence [e.g., Matthaeus and Lamkin, 
1986; Lazarian and Vishniac, 1999; Loureiro et al., 2009; Hi¬ 
gashimori et al., 2013; Yokoi et al, 2013]. In addition to its 
role in the dynamic evolution of reconnection, turbulence is 
believed to play an essential role not only in thermal plasma 
heating but also in nonthermal particle production through 
the stochastic scattering of particles [e.g., Veltri et al., 1998; 
Greco et al, 2002; Zelenyi et al, 1998, 2011; Lazarian et al., 
2012]. However, in previous studies, it was postulated that 
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turbulence can be generated in the high /3 plasma sheet with 
a high magnetic Reynolds number, but the detailed mech¬ 
anism required to generate Alfvenic waves and turbulence 
during the dynamic evolution of reconnection is not under¬ 
stood as yet. 

Several possible candidates for generating Alfvenic fluc¬ 
tuations in the plasma sheet have been proposed: (1) 
the anisotropic ion beams observed in the plasma sheet 
boundary layer (PSBL) can generate magnetohydrodynamic 
(MHD) waves through a family of ion-ion beam instabilities 
[e.g., Gary, 1991; Krauss-Varban and Omidi, 1995; Grig- 
orenko et al., 2011]. The excited Alfvenic wave in the bound¬ 
ary may penetrate into the plasma sheet because the refrac¬ 
tive index of the Alfven wave in the plasma sheet is higher 
than in the lobe plasma region. (2) An elongated magnetic 
diffusion region with strong electric currents may make the 
tearing mode unstable, and many small-scale plasmoids cre¬ 
ated by the tearing mode instability can emanate from the 
diffusion region [e.g., Loureiro et al., 2007; Samtaney et al., 
2009; Bhattacharjee et al., 2009; Pucci and Velli, 2014]. (3) 
The magnetic field pile-up region, where the reconnection 
jet stops suddenly, is expected to excite various plasma 
waves by releasing the bulk flow energy [e.g., Hoshino et al., 
1998, 2001]. In dawn-dusk direction parallel to the electric 
current, the interchange instability can be excited, because 
the gradient of the plasma density is opposite to the de¬ 
celerating plasma flow profile, [e.g., Nakamura et al., 2002; 
Lapenta and Bettarini, 2011]. (4) Turbulence in the solar 
wind may penetrate into the magnetotail across the magne¬ 
topause. Even if the magnetosphere is encircled by a closed 
magnetic field line, fluctuations in the magnetosonic waves 
can propagate perpendicular to the magnetic held. (5) In 
addition to the above possible origins of Alfvenic fluctua¬ 
tions, a reconnection jet at Alfvenic speed can be a free 
energy source for the MHD instability from the release of 
bulk flow energy. 

In this paper, we focus on an MHD instability triggered 
by a bulk How plasma in the reconnection downstream, and 
propagating to the parallel to the bulk flow. In this situa¬ 
tion, some reconnection simulations using the hybrid codes 


1 



X - 2 


MASAHIRO HOSHINO: MHD INSTABILITY IN RECONNECTION JETS 


suggested the generation of a kink-type motion of the re¬ 
connection exhaust as one of possible origins of turbulence 
in the plasma sheet [Lottermoser et al., 1998; Arzner and 
Scholer, 2001; Higashimori and Hoshino, 2012; Liu et al., 
2012], and a tearing mode MHD simulation under the local¬ 
ized bulk flow claimed the rapid excitation of the sausage- 
type fluctuations [Sato and Walker, 1982]. In addition to 
these nonlinear simulation results, the MHD linear insta¬ 
bility with the localized bulk flow oriented parallel to the 
neutral sheet, which focuses on only the reconnection ex¬ 
haust, has been extensively investigated by Shigeta et al. 

1985]; Wang et al. [1988]; Lee et al. [1988]; Biskamp et al. 

1998]. 

In the MHD linear stability analysis, three important 
modes have been discussed. One of the unstable modes is the 
“streaming tearing” mode [Shigeta et al, 1985; Wang et al, 
1988], categorized as a resistive MHD instability, and a sec¬ 
ond mode is the “streaming sausage” mode [Lee et al., 1988; 
Biskamp et al, 1998], categorized as an ideal MHD insta¬ 
bility. The eigenfunctions of the streaming tearing/sausage 
modes are similar to the standard tearing mode [Furth et al., 
1963], whose density perturbation is symmetrical to the neu¬ 
tral sheet. The streaming tearing/sausage modes are not 
zero-frequency modes, and the tearing islands drift with the 
bulk flow speed. In addition to the two symmetric unsta¬ 
ble modes above, an asymmetric unstable mode can be ob¬ 
served, which we call the “streaming kink” mode [Lee et al., 
1988; Biskamp et al., 1998]. This mode can be categorized as 
an ideal MHD instability, and has an instability mechanism 
that is similar to the streaming sausage mode. (Note that 
the terminology of the drift-kink instability whose mode is 
propagating along the electric current direction is different 
from the streaming kink instability discussed in this paper 
[Pritchett et al., 1996; Zenitani and Hoshino, 2005; Fuji- 
moto, 2011 ].) 

In this paper, using the standard method of matrix eigen¬ 
value analysis, we studied the above three different types 
of streaming modes from sub-Alfvenic to super-Alfvenic 
regime, by paying special attention to the earth’s magne¬ 
totail where the lobe density is low, because the dependence 
of bulk flow speed on the linear growth rate has not been sys¬ 
tematically investigated, and because the dilute lobe plasma 
environment has not been taken into account. We observed 
that the streaming sausage and streaming kink modes be¬ 
came unstable in the magnetotail when the bulk flow speed 
reaches a value of several tens of percent of the Alfven speed. 
More importantly, we studied the streaming modes under 
the effects of the Hall current and the guide magnetic field, 
and propose a possible origin of turbulence in the plasma 
sheet during magnetic reconnection. In Section 2, the basic 
model is described with a bulk flow plasma in the center 
of the plasma sheet, and we introduce the method used to 
solve the linear stability of the streaming plasmas. In Sec¬ 
tion 3, we discuss our linear stability results for both the 
streaming tearing/sausage modes and the streaming kink 
mode. In Section 4, we summarize our results and provide 
a perspective on the streaming instabilities. 

2. Streaming Plasma Sheet Model and 
Linear Analysis 


the Harris solution, where the magnetic field B x is given by 
B x = Biobe tanh(j//A), (1) 


and By = 0 and B z = const., where A is the thickness of the 
plasma sheet. The plasma density, p, is expressed by 


P = 


Pps _ 

cosh 2 (y/ A) 


Plobe* 


( 2 ) 


In this paper, we assumed a finite background density, piobe, 
exists over the entire plasma sheet. The temperature for the 
uniform component, piobe, was assumed to be Ti 0 b e = 0 . 

In addition to the standard Harris state, we added a finite 
bulk flow, whose velocity profile, v x , is given by 


_ Tjet _ 

COSh 2 (j//Ajet) ’ 


(3) 


and v y = v z = 0 , where Aj e t is the thickness of the bulk 
flow. The bulk flow may mimic the reconnection exhaust 
with an Alfvenic jet. An MHD instability coupled with the 
streaming plasma may be activated because of this stream¬ 
ing plasma. 

We used a set of compressible MHD equations 


dp 

dt 


-V • (pv), 


(4) 


^ + v -W = --Vp+ — jx B + uV 2 v, (5) 
ot p pc 

®=Vx(iixB-i/xB) + ??V 2 B, (6) 

where the terms e and n in the j x B Hall component repre¬ 
sent the electric charge and the number density, respectively. 
The equation of state is assumed to be adiabatic, i.e., p oc p 7 , 
where the ratio of the specific heat is 7 = 5/3. The electron 
temperature was set to zero. The terms v and r/ are the 
fluid viscosity and the electric resistivity, respectively. 

After linearizing the above equations using the standard 
perturbation method in x and y two-dimensional space, 
we solved the set of linearized equations using the matrix 
eigenvalue method, which is often used in the analysis of 
plasma instabilities [e.g., Hoshino, 1991]. The time deriva¬ 
tive, df(x,y,t)/dt, of the physical quantity, f, and the spa¬ 
tial derivative along the plasma sheet, df(x, y, t)/dx, can be 
expressed by —itof(k, y, u>) and ikf(k, y, to) using a Fourier- 
Laplace transformation in time and space, respectively. In 
this paper, we assumed that the physical quantities were 
uniform in the z direction, i.e., d/dz = 0 . 

The spatial derivatives in the direction vertical to the 
plasma sheet, df /dy, can be approximated by the fourth- 
order, finite difference representation of the differentiation 
with respect to y. In addition to this, we used a nonuniform 
grid spacing. For example, the first derivative in the finite 
difference representation df /dy can be expressed by 


In this paper, we discuss the streaming MHD modes in the 
framework of the Hall MHD. As discussed later, the stream¬ 
ing instability itself can be activated in the standard MHD 
regime without the Hall effect, but we use the Hall MHD 
equation having its application to the earth’s magnetotail 
in mind, because from satellite observations, the thickness 
of the reconnection region is known to be of the order of the 
ion inertial length [e.g., Sergeev et al., 1993; Asano et al., 
2003; Nakamura, 2006]. 

We studied the linear stability of a Harris-type plasma 
sheet with a finite bulk flow oriented parallel to an antipar- 
allcl magnetic field. The initial equilibrium state was set to 
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_ Aj_i + A(_2 Aj+i + Ai +2 
Ai_iA,_2 Ai+iAi +2 ’ 

where A 1± i = \yi±i - yi\ and A i±2 = \yt ±2 - Vi\- A nonuni¬ 
form grid size is useful for resolving the resistive layer in the 
vicinity of the neutral sheet for the case of a large magnetic 
Reynolds number [e.g., Hoshino, 1991]. The grid size in the 
vicinity of the neutral sheet was set to be five to 50 times 
smaller than that in the lobe region. The number of grid 
points was chosen to be in the range 400 to 1600 to ensure 
the convergence of the numerical results. 

The size of the plasma sheet, \y\ < L , was set to be 
L/A = 5, and the boundary condition was assumed to 
be f(x, \y\ = L,t) = 0 for the spontaneous reconnection 
model. In our study, we used the magnetic Reynolds num¬ 
ber, Rm = VaX/t/ = 10 3 , and the fluid Reynolds number, 
Ry = VaX/v = 10 4 , in most of our calculations, except 
for the cases mentioned. The other plasma parameters dis¬ 
cussed in this paper are listed in Table 1. 

3. Results of the Streaming MHD Instability 

3.1. Streaming instability with symmetric or 
asymmetric perturbation 

Figures la-c show the linear growth rates obtained from 
our matrix eigenvalue analysis for three different bulk flow 
speeds. The top of each panel shows the linear growth rate, 
Im(u;TA), while the bottom of each panel shows the oscil¬ 
lation frequency, Re(wTA), where ta = A /Va is the Alfven 
transit time. Figure la shows the case where there was no 
bulk flow, Vjet = 0 , and we only obtained the standard tear¬ 
ing mode [Furth et al., 1963], where the unstable region only 
appeared for fcA < 1, and the oscillation frequency was zero. 

Figure lb shows the case for a finite bulk flow speed with 
Vjet = fjet /Va, lobe = 0 . 2 , normalized to the lobe Alfven 
speed, defined by VA.iobe = 5ic.be/V^Piobe- Note that the 
bulk flow speed is normalized to the lobe Alfven speed, while 
the Alfven transit time was normalized to the Alfven speed, 
defined by the total plasma density in the neutral sheet, 
po = Pps + piobe, because the reconnection jet speed was ex¬ 
pected to reach to the lobe Alfven speed when a switch-off 
shock formed in the boundary. 

In Figure lb, it can be seen that the unstable region ex¬ 
tended to the larger wave number region for fcA > 1, and 
that the growth rate was enhanced slightly compared with 
the standard tearing mode in Figure la. The oscillation 
frequency, Re(wTA), denoted by the red-colored circles can 
be approximated by the drift frequency/Doppler shift fre¬ 
quency of Re(wrA) — (kA)Vjet, depicted by the black curve. 
We confirmed that the behavior of the linear growth rate was 
basically the same as the “streaming tearing” mode studied 
by Shigeta et al. [1985]; Wang et al. [1988], and we found 
that from the eigenfunctions the density perturbation was 
symmetric to the neutral sheet (not shown here). 

Figure lc shows the growth rate and the oscillation 
frequency for a higher bulk flow speed with Vj e t = 
Ujet/VA.iobe = 0.4. In addition to the streaming tearing 
mode denoted by the red-colored circles, we observed the 
streaming kink mode, denoted by the green-colored squares, 
whose density perturbation, p, and vector potential, A z , be¬ 
came asymmetric versus the neutral sheet. The streaming 
kink mode could not be excited below a given minimum 
threshold speed. However, since the threshold speed was 
less than the lobe Alfven velocity, VA.iobe, the streaming 
kink mode can be easily generated during magnetic recon¬ 
nection. Another interesting point is that the growth rate 
of the streaming kink mode was larger than the streaming 
tearing mode, and the streaming kink mode may dominate 
the plasma sheet turbulence in the reconnection jet. 

Figure 2 shows typical linear structures of the vector 
potential, A~, and plasma density, p, in (x, y) space that 
were reconstructed from the eigenvalues and eigenfunctions 


in Fourier ( k, y) space. The first-order perturbations were 
superposed onto the zero-order quantities, and we assumed 
the magnitude of the first-order peak density was 30% of the 
zero-order density in the neutral sheet. Figure 2a shows the 
streaming tearing mode for (Vj e t,fcA) = (0.2, 0.3); Figure 2 
b shows the streaming kink mode for (Vj e t, kX) = (0.4,1.5). 
The structure of the streaming tearing mode was similar to 
that of the standard tearing mode, except for the magnetic 
islands that streamed with the zero-order bulk flow speed, 
Ujet- On the other hand, the structure of the streaming kink 
mode had a periodic distortion in the neutral sheet. 

3.2. Maximum growth rate 

So far, we have discussed that the streaming tearing mode 
is unstable for a relatively slow bulk flow regime, while the 
streaming kink mode appears to be unstable in a relatively 
fast flow regime. Let us now study the behavior of the two 
unstable modes as a function of the bulk flow speed, Vj e t- 

Figure 3 shows the maximum growth rate and the corre¬ 
sponding wave number as a function of the bulk flow speed, 
Vjet - In our numerical calculations, we surveyed the growth 
rate and the oscillation frequency for a given wave number 
for fcA = 0.05 x i , for i = 0 to 63. The wave number and bulk 
flow were normalized to the thickness of the plasma sheet, A, 
and the lobe Alfven velocity, VA.iobe, respectively. The red- 
colored circles denote the streaming tearing mode, while the 
green-colored circles denote the streaming kink mode. We 
observed that the streaming tearing mode with a symmetric 
density perturbation dominated the regime where Vj e t < 0.4, 
while the growth rate of the streaming kink mode became 
larger than that of the streaming tearing mode when the 
bulk flow speed exceeded about Vjet > 0.4. These modes 
have a long wavelength nature for fcA < 2. 

By carefully examining the streaming tearing mode de¬ 
noted by the red-colored circles, we could observe that 
the characteristics of the unstable wavelengths were differ¬ 
ent below and above a value of Vj e t ~ 0.25. The unsta¬ 
ble wavelengths for Vjet < 0.25 were the long wavelengths 
with fcA < 1, while for Vj e t > 0.25, the unstable wave¬ 
lengths occurred for 1 < fcA < 2. Moreover, the growth 
curve of Im(u;TA) showed a discontinuous transition from 
Im(urrA) ~ 10~ 2 to 5 x 10~ 2 around a value of Vj e t ~ 0.25. 
This suggests that two different modes may exist below and 
above V]et ~ 0.25. 

As we discuss this behavior further in the next subsection, 

3.3. it can be approximated that the linear growth rate of the 
symmetric unstable mode for Vj e t < 0.25 strongly depends 
on the magnetic Reynolds number, while the linear growth 
rate of the symmetric unstable mode for Vj e t > 0.25 does 
not depend on the magnetic Reynolds number. Therefore, 
we distinguish the symmetric perturbations from their de¬ 
pendence of the linear growth rate on the magnetic Reynolds 
number, Rm- One of these modes is the “streaming tear¬ 
ing” mode, categorized as a resistive MHD mode [Shigeta 
et al. , 1985; Wang et al, 1988], and the other mode is the 
“streaming sausage” mode, classified as an ideal MHD mode 
[Lee et al ., 1988; Biskamp et al., 1998]. The streaming tear¬ 
ing and streaming sausage modes are observed in relatively 
slow and fast bulk flow speeds, respectively. 

The perturbed structure of the streaming sausage mode 
is similar to that of the streaming tearing mode (not shown 
here), and the vector potential, A z , and the plasma den¬ 
sity, p, for the streaming tearing/sausage modes show sim¬ 
ilar properties to the standard tearing mode, namely those 
perturbations are symmetric against the neutral sheet. 

3.3. Reynolds number dependence, Rm 
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Figure 4 shows the linear growth rates for the stream¬ 
ing tearing, streaming sausage, and streaming kink modes 
as a function of the magnetic Reynolds number. The bulk 
flow speed was fixed to be Vj e t =0.5 for all cases. The red- 
colored circles with the solid line and the red-colored circles 
with the dashed line denote the growth rates of the symmet¬ 
ric density perturbation for the streaming tearing mode with 
kX = 0.3 and the streaming sausage mode with kX = 1.0, 
respectively. The green-colored squares correspond to the 
streaming kink mode with an asymmetric density perturba¬ 
tion for the wave number with kX = 1. There may be more 
than one unstable mode for any given wave number. The 
streaming sausage and streaming kink modes showed sym¬ 
metric and asymmetric density fluctuations, respectively. 

The growth rate of the streaming tearing mode decreased 
with increasing magnetic Reynolds number, Rm, and the 
growth rate was approximately Im(a;rA) oc R~^^ 3 , which 
showed a similar dependence to that of the standard tear¬ 
ing mode of R^ 3 3 [Furth et al, 1963], in the sense that 
it was a decreasing function. The streaming tearing mode 
is slightly enhanced for a large magnetic Reynolds number, 
by utilizing the free energy of the streaming bulk flow. On 
the other hand, the growth rate of the streaming sausage 
and streaming kink modes did not depend strongly on the 
Reynolds number, and the growth rates were almost con¬ 
stant in the high magnetic Reynolds number regime. This 
result is consistent with the streaming sausage and stream¬ 
ing kink modes under an ideal MHD regime discussed by Lee 
et al. [1988]; Bishamp et al. [1998]. Therefore, in a plasma 
medium with a high magnetic Reynolds number, Rm, the 
streaming sausage and streaming kink modes may play an 
important role in the excitation of Alfvenic fluctuations, 
even if the streaming tearing mode cannot be excited. 

3.4. Lobe density dependence, piobe/Po 

We confirmed that the streaming tearing/sausage and 
streaming kink instabilities can be excited by adding a bulk 
flow plasma to the plasma sheet [e.g., Shigeta et al., 1985; 
Wang et al., 1988; Lee et al., 1988; Biskamp et al., 1998], 
and these instabilities are thought to be important in the 
reconnection jet region. Since the reconnection jet speed 
is known to be the lobe Alfven velocity if a switch-off slow 
mode shock is formed, the generation of turbulence by the 
streaming instabilities may be strongly controlled by the 
lobe plasma density. 

To study the dependence of the lobe plasma density on 
the linear stability, the data in Figure 5 show the growth rate 
for piobe/po = 0.01, 0.05, and 0.1, where po = p ps + piobe is 
the total plasma density in the neutral sheet [e.g., Ishisaka 
et al., 2001]. Figure 5a shows the same data as in Figure 
3. It can be seen that the streaming tearing and streaming 
sausage modes with a symmetric perturbation are excited 
in the relatively slow bulk flow regime, while the streaming 
kink mode with an asymmetric perturbation appears in a 
relatively fast flow regime. The onset speeds of the stream¬ 
ing sausage and kink modes increased with increasing lobe 
plasma density, piobe/po- 

Discontinuous features in the growth rates of the symmet¬ 
ric perturbation, suggesting a transition from the stream¬ 
ing tearing mode to the streaming sausage mode, were ob¬ 
served around Vj e t ~ 0.25, 0.55, and 0.75 for piobe/po = 0.01, 
0.05, and 0.1, respectively. The onset speed of the transi¬ 
tion also increased with increasing the lobe plasma density. 
The growth rates for the streaming tearing and streaming 
sausage modes did not change appreciably on increasing the 
lobe plasma density, but growth rates for the streaming kink 
mode could be suppressed because of the contribution of the 
dense lobe plasma. 

Let us now look at the detail on the streaming tearing 
mode. We can recognize modest peaks in the growth rates 
for relatively small jet velocities, namely the positions of 


the peaks are Vj e t ~ 0.1 for piobe/po = 0.01, Vj e t ~ 0.2 for 
piobe/Po = 0.05, and V iet ~ 0.3 for pi obe /po = 0.1. We find 
these peaks almost co rrespond to the sound speeds in the 
plasma sheet, c s = y^TpsPps/po, which are c s /V A , lobe = 
0.091, 0.204, and 0.289 for pi ob e/po = 0.01, 0.05, and 0.1, 
respectively. The streaming tearing mode can be amplified 
in the subsonic jet flow, while it can be suppressed in the 
supersonic jet flow. 

The reason why the growth rate of the streaming tear¬ 
ing instability is enhanced in the subsonic flow may be in¬ 
terpreted from the decrease in gas pressure around the X- 
type point, and the resulting emission of a fast expansion 
wave. The streaming plasma around the X-type point can 
be squeezed by the reconnecting magnetic field lines, and as 
a result, the streaming plasma speed can increase. Based 
on Bernoulli’s principle of the dynamic and gas pressure 
balance along a streamline with pv 2 / 2 + p = const., the gas 
pressure decreases, and the inflow plasma toward the X-type 
region can be enhanced. On the other hand, for the case of 
a supersonic flow, the streaming plasma flowing around the 
squeezed X-type point can be decelerated, and then the gas 
pressure increases. Therefore, the inflow of plasma toward 
the X-type point can be suppressed. 

A transition from streaming tearing to streaming sausage 
modes was observed around Vj e t = 0.25, 0.55, and 0.75 for 
piobe /Po = 0.01, 0.05, and 0.1, respectively. The streaming 
sausage and streaming kink modes are classified as nonresis- 
tive MHD instabilities, and the mechanism of the instability 
may be controlled by the balance between the tension force 
of the magnetic field and the centrifugal force of the bulk 
flow under an infinitesimal distorted magnetic field line. If 
the bulk flow speed is fast enough to overcome the tension 
force of the magnetic field line, then we expect the stream¬ 
ing sausage and streaming kink instabilities to occur. Let us 
assume that the distortion of the magnetic field lines occurs 
at the boundary between the lobe and the plasma sheet at 
y = aAjet, where a ~ 1 is a tuning parameter determined by 
the details of the plasma process. Then, the marginal state 
of the instability can be expressed by 


47T 


\ cosh 2 (o') 


where a' = aXj et /X. In this case, the onset velocity of the 
streaming sausage may be given by 


Vjet =cosh 2 (a)tanh(a')/wU- \ ■ (9) 

V COSll(a ') 2 piobe 

For a = a' = 0.96 with Aj e t = A, we obtain Vj e t = 0.25, 0.54, 
and 0.75 for piobe/po = 0.01, 0.05, and 0.1, respectively. 
This simple estimation can model all the transition veloc¬ 
ities from the streaming tearing to the streaming sausage 
modes very well. 

The transition from the streaming sausage mode to the 
streaming kink mode, which can be recognized by the change 
in symbol from the red-colored circles to the green-colored 
squares in Figure 5, appears around Vj e t = 0.4, 0.8, and 
1.05 for piobe/Po = 0.01, 0.05, and 0.1, respectively. These 
transition velocities can be modeled by assuming a = 1.1 in 
Eq. (9), and again, we obtain Vj e t = 0.36, 0.80, and 1.08 for 
piobe/po = 0.01, 0.05, and 0.1, respectively. 

In the earth’s magnetotail, the lobe plasma density is 
lower than the plasma sheet density, i.e., piobe/p ps ~ 0.01 — 
0.05 [e.g. Ishisaka et al., 2001]. On the other hand, many 
MHD simulation studies of magnetic reconnection assume 
a relatively higher plasma density exists in the lobe to re¬ 
duce the Alfven velocity and to save on computational CPU 
time. A high density plasma can suppress the emission of 
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the streaming kink instability, and this may be one reason 
why no streaming kink instability in association with MHD 
magnetic reconnection simulations has been reported. 

3.5. Effect of bulk flow size, Aj e t 

It is also of interest to study the effect of the jet size Aj e t, 
because the free parameter Aj e t is independent of the equi¬ 
librium state. Figure 5d shows the result where Aj e t = 0.5, 
but keeping piobe /po = 0.1, the same as in Figure 5c. By 
comparing Figures 5c and 5d, we find that the growth rate 
of the streaming tearing mode can be enhanced in the range 
Vjet < 0.25, namely almost in the subsonic flow regime, and 
is reduced in the supersonic flow regime, where the sound 
speed is c s /VA,iobe = 0.289. The enhancement and suppres¬ 
sion of the growth rate is much magnified in the case of 
Figure 5d. The reason for this is probably because of the 
localization of the instability inside the high /? plasma sheet, 
where the magnetic field is weak. That is, for the case of a 
narrow jet, the instability appears close to the neutral sheet 
with a weak magnetic field, and the effect of the magnetic 
tension/pressure becomes less important compared with the 
plasma dynamic and gas pressures. 

The onset velocity of the streaming sausage and stream¬ 
ing kink modes for a narrow jet is reduced compared with 
Figure 5c. This can also be interpreted by the weak mag¬ 
netic field tension force. By using the modeling of the onset 
jet velocities in Eq. (9) with the same values of a = 0.96 and 
a' = 0.96(Aj e t/A) = 0.48, we obtained Vj e t = 0.45, which 
agrees very well with the onset of the streaming sausage 
mode. For the transition from the streaming sausage mode 
to the streaming kink mode, by substituting the same value 
of a = 1.1 and a' = l.l(Aj e t/A) = 0.55, we obtained 
Vjet = 0.68. Again, we obtain a good agreement between 
the linear analysis and the theoretical modeling. 

The enhancement of the linear growth rate can be under¬ 
stood from the localization of the instability in the vicinity 
of the neutral sheet because of the narrow jet flow. For 
a distorted magnetic field line, the magnetic tension force 
becomes weak, while the centrifugal force of the bulk flow 
remains constant, if the bulk flow speed is constant. There¬ 
fore, the linear growth rate can be enhanced by reducing 
the size of the bulk flow. However, since the total free en¬ 
ergy required to excite the streaming instabilities is reduced, 
the magnitude of the amplitude of the nonlinear saturation 
of the streaming instabilities may remain small. A discus¬ 
sion of the saturation level is beyond the scope of our linear 
analysis. 

3.6. Mode structure under the Hall effect 

So far, we have discussed the behavior of the streaming 
modes by neglecting the Hall effect, that is, we have as¬ 
sumed that the ion inertia length is much smaller than the 
thickness of the plasma sheet. However, it is known that 
the thickness of the plasma sheet is of the order of the ion 
inertia length near the reconnection region. In this subsec¬ 
tion, we discuss the effect of the Hall term on the streaming 
instabilities. 

Figure 6 shows the maximum growth rate and its corre¬ 
sponding wave number as a function of the bulk flow speed, 
Vjet, under the Hall MHD with (VA/fL)/A = 1.0. The other 
parameters are the same as those shown in Figure 3 with¬ 
out any Hall effect. By comparing Figures 3 and 6, we find 
that: (1) the linear growth rate of the Hall tearing mode 
is larger than that of the standard MHD tearing mode at 
Vjet = 0; (2) regardless of the Hall term, the three different 
streaming tearing, streaming sausage, and streaming kink 
modes are separated by bulk flow speeds of Vj e t ~ 0.2, and 
0.3; and (3) under a Hall effect, the growth rate does not 
change much but the unstable region is shifted to slightly 
longer wavelengths. 

Let us now look at the eigenstructure of the streaming 
modes under the Hall effect shown in Figure 6. Figures 7 a- 
c show the reconstructed structures for Vj e t = 0, 0.2, and 


0.4 in the two-dimensional x-y space, respectively. We have 
chosen the wave numbers with the maximum growth rate 
for the fixed bulk jet speed, namely k\ = 0.25, 0.5, and 0.5 
for Vjet = 0, 0 . 2 , and 0.4, respectively. The reconstructed 
structures were obtained by superposing the first-order per¬ 
turbations onto the zero-order quantities, and we assumed 
the magnitude of the first-order peak density was 30% of the 
zero-order peak density in the neutral sheet. 

In Figure 7 (a), the quadrupole magnetic field structure 
B z from the Hall effect can be seen to be localized in the 
vicinity of the neutral sheet. Note that we assumed a mag¬ 
netic Reynolds number of Rm = 10 3 , but a lower value of 
Rm can produce a wider quadrupole magnetic field struc¬ 
ture. By increasing the bulk flow speed up to (b) Vj e t = 0.2 
and (c) Vjet = 0.4, in addition to the localized B z component 
in the vicinity of the neutral sheet, large-scale fluctuations 
outside the localized B z component for both the stream¬ 
ing sausage mode and the streaming kink mode were clearly 
seen. The polarity of B z is asymmetric for the streaming 
sausage mode and symmetric for the streaming kink mode. 
In association with the generation of B z , the global plasma 
flow pattern in the entire plasma sheet can be seen, as rec¬ 
ognized from the white arrows in the vertical direction. We 
believe that the Hall term may have a significant effect on 
the generation of Alfvenic fluctuations in the plasma sheet. 

3.7. Guide magnetic field effect 

It is also interesting to study the effect of the guide mag¬ 
netic field B z on the streaming instabilities, because a global 
dynamical behavior of magnetic reconnection is known to 
be controlled by the guide magnetic field. In fact, the 
guide magnetic field in the magnetotail, which is the dawn- 
dusk magnetic field, is known to exist due to the penetra¬ 
tion of the solar wind magnetic field [e.g., Cowley , 1981; 
Petrukovich, 2011; Rong et al., 2012]. 

Shown in Figure 8 is the growth rates of the standard 
tearing mode and the streaming tearing, sausage and kink 
modes as a function of the initial guide magnetic field B z . 
The guide magnetic field is normalized by the lobe mag¬ 
netic field .Biobe, and the growth rate is normalized by the 
Alfven transit time. The growth rate of the standard tear¬ 
ing mode with (kX, Vjet) = (0.3,0) denoted by the blue- 
solid line does not depend on the magnitude of the guide 
magnetic field, and that of the streaming tearing mode 
with (fcA,Vj e t) = (0.3, 0 . 8 ) denoted by the blue-dashed line 
slightly increases as increasing the guide magnetic field. 
However, the streaming sausage and kink modes with the 
shorter wavelengths show significant change on the guide 
magnetic field. The red-solid and red-dashed lines show 
respectively the sausage mode with (fcA,Vj e t) = ( 1 . 0 , 0 . 8 ) 
and (/cA, Vjet) = (1.5, 0 . 8 ), while the green-solid and green- 
dashed lines are the kink mode with (fcA,Vj e t) = ( 1 . 0 , 0 . 8 ) 
and (kX, Vjet) = (1.5, 0.8), respectively. These sausage and 
kink modes have a tendency to be stabilized for a relatively 
weak guide magnetic field B z /B\ 0 b e < 0.5, while they are 
destabilized for a large guide magnetic field B z /B\ 0 b e > 0.5. 

It should be noted that the inclusion of the guide mag¬ 
netic field may reduce the plasma compressibility, and that 
the velocity shear instability such as Kelvin-Helmholtz insta¬ 
bility has been shown to be destabilized by a strong vertical 
magnetic field perpendicular to the shear plasma flow [Miura 
and Pritchett, 1982]. The similarity of the basic behav¬ 
ior between Kelvin-Helmholtz instability and the streaming 
sausage and kink instability has been discussed by Lee et al. 
[1988]; Biskamp et al. [1998], and the incompressible limit 
due to the strong magnetic field may enhance the growth 
rates of the streaming sausage and kink instabilities. 
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4. Discussion and Conclusions 

We have observed that the streaming MHD instability can 
be activated by not only a super-Alfenic bulk flow but also 
a sub-Alfenic flow, oriented parallel to the antiparallel mag¬ 
netic field by taking into account of the dilute plasma envi¬ 
ronment in the lobe, and found that the growth rates were 
enhanced more strongly than the standard tearing mode 
without bulk flow. We also studied that the guide mag¬ 
netic field perpendicular to the bulk flow can significantly 
modify the growth rates of the streaming sausage and kink 
modes. Furthermore, we also found that Alfvenic pertur¬ 
bations of B z could be generated outside the plasma sheet 
when the Hall effect is included. The perturbations may 
play an important role in the emission of MHD waves over 
the plasma sheet. We believe that the streaming modes may 
be a candidate for plasma sheet turbulence. 

Since we observed that the bulk flow embedded in the 
plasma flow could excite three different streaming modes, 
depending on the bulk flow speed, the next important is¬ 
sue to understand is the nonlinear evolution and the satu¬ 
ration of the instabilities. Sato and Walker [1982] investi¬ 
gated the streaming tearing mode using MHD simulations, 
and suggested that the streaming tearing mode is excited 
much more in the presence of bulk flow. However, as Sato 
and Walker [1982] assumed that the initial plasma density 
was uniform over the plasma sheet, they only discussed the 
streaming tearing mode with a symmetric density perturba¬ 
tion against the neutral sheet. Because a low density plasma 
in the lobe region is required to excite the kink mode, based 
on our linear stability analysis, a streaming kink mode was 
not demonstrated in their simulations. It is important to in¬ 
vestigate the nonlinear evolution of the streaming kink mode 
as well. Interestingly, some hybrid simulations with a low 
density plasma in the lobe have suggested the presence of a 
streaming kink-like structure in the reconnection jet [Lotter- 
moser et al., 1998; Arzner and Scholer, 2001; Higashimori 
and Hoshino, 2012; Liu et al., 2012]. However, since the lobe 
densities in these simulations are still higher than that in the 
earth’s magnetotail, we expect that a streaming kink-like 
structure can be excited easily in more realistic situations. 

In this paper, we have only discussed the linear stability 
of the streaming MHD instability for an isotropic pressure, 
i.e., P|| = p± ; however, under a collisionless reconnection, 
it is known that the anisotropy in pressure for py > p± 
can be generated [e.g., Hirabayashi and Hoshino, 2013]. In 
this situation, the fire-hose instability may couple with the 
streaming modes [ Arzner and Scholer, 2001]. Under the 
anisotropic plasma situation with p\\ > p±, the magnetic 
tension force for the distorted magnetic field line becomes 
weak; as a result, we may expect that the streaming instabil¬ 
ity can be easily excited. However, a study of the behavior 
of an anisotropic plasma is beyond the scope of this paper. 

We assumed that the physical quantities were uniform in 
the z direction, i.e., d/dz = 0, and that the Fourier modes 
are restricted to be propagating parallel to the x direction. 
However, it would be interesting to understand the oblique 
propagation as well, because the drifting plasma population 
parallel to the electric current may couple with the stream¬ 
ing plasma flow. In addition, if there is the density gradient 
along the plasma flow, the interchange instability may be ex¬ 
cited in the z direction [e.g., Nakamura et al., 2002; Lapenta 
and Bettarini, 2011]. The three dimensional behavior of the 
streaming instability would play an important role on the 
actual plasma sheet. 

We have discussed the streaming tearing, streaming 
sausage, and streaming kink modes as possible models of 
generation of MHD turbulence in the earth’s magnetotail. 
However, the streaming instability may be found in other 
applications, one of which is the so-called “channel flow” 


formed during the magnetorotational instability (MRI) for 
accretion disks [e.g., Balbus and Hawley, 1991]. It is of inter¬ 
est to note that the particle-in-cell simulation results on the 
MRI by Hoshino [2013] show a streaming kink-type channel 
flow occurring just before the onset of reconnection. The 
streaming kink mode may also play an important role on 
the mass and angular momentum transport during the for¬ 
mation of an accretion disk. 
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Table 1 . Plasma parameters used in our linear analysis: the 
normalized lobe plasma density N\ 0 ^ e = pi o b e /P0i the nor¬ 
malized bulk flow speed by the lobe Alfven speed Vj e t, the 
normalized wavelength k\, the ratio of the bulk flow jet size 
Aj et and the thickness of the plasma sheet A, the Hall effect 
on (VA/Oi)/A, and the guide magnetic field TL/Biobe- 
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Figure 1. Linear growth rate, Im(tUTA), (top) and oscil¬ 
lation Frequency, Re(cPTA), (bottom) as the function of 
wave number, kX, for three different bulk flow speeds: (a) 
Vjet = Pjet/VA,lobe = 0, i.e., nO bulk flow, (b) Vjet = 0.2, 
and (c) Vjet — 0.4. The solid black lines in the bottom 
panels (b) and (c) are the Doppler shift frequencies of 
Re(urrA) = (fcA)Vjet. 
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Figure 2. Unstable structures for the vector potential, 
A z , (top) and the plasma density, p, (bottom), which 
were reconstructed from the eigenfunctions and eigenval¬ 
ues. The flow vectors (red-colored arrows) are super¬ 
posed on the magnetic field lines (black) of the contours 
for A z . (a) The streaming tearing mode with Vj e t = 0.2 
and kX = 0.3 in Figure lb, and (b) the streaming kink 
mode with V) e t = 0.4 and kX = 1.5 in Figure lc. 
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Figure 3. Maximum growth rate (top) and the corre¬ 
sponding wave number (bottom) as a function of the bulk 
flow speed, Vj e t, for (Va/^i) / X = 0 without a Hall effect. 
The red-colored circles and the green-colored squares 
show the symmetric and asymmetric density perturba¬ 
tions, respectively. 
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Figure 4. Growth rates, Im(a>TA), for the streaming 
tearing mode for fcA = 0.3 (red-colored solid line with 
circles), the streaming sausage mode for fcA = 1.0 (red- 
colored dashed line with circles), and the streaming kink 
mode for fcA = 1.0 (green-colored solid line with squares) 
as a function of magnetic Reynolds number, Rm- The 
bulk flow speed was set to be Vj e t = 0.5 for all modes. 
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Figure 5. Dependence of the lobe plasma density, 
Piobe/Pps, and the bulk flow size, Aj e t. Maximum growth 
rates, Im(u;TA), are shown as a function of the bulk flow 
speed, Vjet, for the following cases: (a) piobe/po = 0.01; 
(b) piobe/po = 0.05; and (c) piobe/po = 0.1. (d) The lobe 
density is the same as in case (c), but the thickness of 
the bulk flow jet is narrower, with Aj e t/A = 0.5. The 
red-colored circles and the green-colored squares show 
the symmetric (tearing/sausage) and asymmetric (kink) 
density perturbations, respectively. 
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Figure 6. Maximum growth rate (top) and its corre¬ 
sponding wave number (bottom) are shown as a func¬ 
tion of the bulk flow speed, under the Hall effect with 
(Va/ 1L)/A = 1. Except for the Hall parameter, the other 
parameters are the same as those used in Figure 3. 



Figure 7. Unstable structures for the plasma density, p, 
(top), the flow vectors, v, (white arrows), and the mag¬ 
netic field B z (bottom): (a) the standard tearing mode 
without a bulk flow speed Vj e t = 0 and kX = 0.3; (b) the 
streaming sausage mode with Vj e t = 0.2 and k\ = 0.5; 
and (c) the streaming kink mode with Vj e t = 0.4 and 
kX = 0.55. 
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Figure 8. Growth rates as a function of the guide 
magnetic field for four different unstable modes: (a) the 
standard tearing mode with (AA, Vj e t) = (0.3,0) (blue 
and solid line), (b) the streaming tearing mode with 
(A:A, Vjet) = (0.3,0.8) (blue and dashed line), (c) the 
streaming sausage mode with (AA,Vj et ) = (1.0, 0.8) (red 
and solid line) and with (AA, Vj e t) = (1.5, 0.8) (red and 
dashed line), and (d) the streaming kink mode with 
(AA, Vjet) = (1.0,0.8) (green and solid line) and with 
(AA, Viet) = (1.5,0.8) (green and dashed line). 





